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Abstract 

In this this paper the stochastic background of gravitational waves 
(SBGWs) is analyzed with the auxilium of the WMAP data. We empha- 
size that, in general, in previous works in the literature about the SBGWs, 
old COBE data were used. After this, we want to face the problem of how 
the SBGWs and f(R) gravity can be related, showing, vice versa, that a 
revealed SBGWs could be a powerly probe for a given theory of gravity. 
In this way, it will also be shown that the conformal treatment of SBGWs 
can be used to parametrize in a natural way f(R) theories. 
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1 Introduction 



The accelerate expansion of the Universe, which is today observed, shows that 
cosmological dynamic is dominated by the so called Dark Energy which gives 
a large negative pressure. This is the standard picture, in which such new in- 
gredient is considered as a source of the rhs of the field equations. It should 
be some form of un-clustered non-zero vacuum energy which, together with the 
clustered Dark Matter, drives the global dynamics. This is the so called "con- 
cordance model" (ACDM) which gives, in agreement with the CMBR, LSS and 
SNela data, a good trapestry of the today observed Universe, but presents sev- 
eral shortcomings as the well known "coincidence" and "cosmological constant" 
problems [JJ. 

An alternative approach is changing the Ihs of the field equations, seeing if 
observed cosmic dynamics can be achieved extending general relativity (2j [3J HI 
\E\. In this different context, it is not required to find out candidates for Dark 
Energy and Dark Matter, that, till now, have not been found, but only the 
"observed" ingredients, which are curvature and baryonic matter, have to be 
taken into account. Considering this point of view, one can think that gravity 
is not scale- invariant [5] and a room for alternative theories is present [3 \7\ [8] . 
In principle, the most popular Dark Energy and Dark Matter models can be 
achieved considering f(R) theories of gravity [21 Gil 0], where R is the Ricci 
curvature scalar. In this picture even the sensitive detectors for gravitational 
waves, like bars and interferometers (i.e. those which are currently in operation 
and the ones which are in a phase of planning and proposal stages) [3 [13 [Til 
[121 [TH [HI HH [16] , could, in principle, be important to confirm or ruling out the 
physical consistency of general relativity or of any other theory of gravitation. 
This is because, in the context of Extended Theories of Gravity, some differences 
between General Relativity and the others theories can be pointed out starting 
by the linearized theory of gravity [3 [17\ HH [19] . 

This philosophy can be taken into account also for the SBGWs which, to- 
gether with cosmic microwave background radiation (CMBR), would carry, if 
detected, a huge amount of information on the early stages of the Universe 
evolution [21 [23 El])- Also in this case, a key role for the production and the 
detection of this graviton background is played by the adopted theory of gravity 

EJj. 

In the second section of this paper the SBGWs is analyzed with the auxilium 
of the WMAP data [231221122]- We emphasize that, in general, in previous works 
in the literature about the SBGWs, old COBE data were used (see [21 [23 (23 
[221 [281 [29] for example). 

In the third section we want to face the problem of how the SBGWs and 
f(R) gravity can be related, showing, vice versa, that a revealed SBGWs could 
be a powerly probe for a given theory of gravity. In this way, it will also be 
shown that the conformal treatment of SBGWs can be used to parametrize in 
a natural way f(R) theories. 
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2 Tuning the stochastic background of gravita- 
tional waves using the WMAP data 

From our analysis, it will result that the WMAP bounds on the energy spectrum 
and on the characteristic amplitude of the SBGWs are greater than the COBE 
ones, but they are also far below frequencies of the earth-based antennas band. 
At the end of this section a lower bound for the integration time of a potential 
detection with advanced LIGO is released and compared with the previous one 
arising from the old COBE data. Even if the new lower bound is minor than 
the previous one, it results very long, thus for a possible detection we hope in 
the LISA interferometer and in a further growth in the sensitivity of advanced 
projects. 

The strongest constraint on the spectrum of the relic SBGWs in the fre- 
quency range of ground based antennas like bars and interferometers, which is 
the range 10Hz < f < l0 4 Hz, comes from the high isotropy observed in the 
CMBR. 

The fluctuation AT of the temperature of CBR from its mean value To = 
2.728 K varies from point to point in the sky (20, 22, 23], and, since the sky can 
be considered the surface of a sphere, the fitting of AT is performed in terms 
of a Laplace series of spherical harmonics 

AT °° 1 

-t-@) = E E ai«3Mft)» (!) 

l — l m——L 

and the fluctations are assumed statistically independent (< ct; TO >= 0, < 
ai m a\, m i >= CiSu'Smm')- In eq. Ct denotes a point on the 2-sphere while the 
ai m are the multipole moments of the CMBR. For details about the definition of 
statistically independent fluctations in the context of the temperature fluctations 
of CMBR see [2g|l25]. 

The WMAP data [22l E3] permit a more precise determination of the rms 
quadrupole moment of the fluctations than the COBE data 

Q rms ee T( V J^)i = 8 ± IpK, (2) 

z — ' 47T 

m=— 2 

while in the COBE data we had [25l [26] [27] 

Q rms = lAJ&tllnK. (3) 

A connection between the fluctuation of the temperature of the CMBR and the 
SBGWs derives from the Sachs-Wolfe effect [12 [30]. Sachs and Wolfe showed 
that variations in the density of cosmological fluid and GWs perturbations re- 
sult in the fluctuation of the temperature of the CMBR, even if the surface 
of last scattering had a perfectly uniform temperature [30]. In particular the 
fluctuation of the temperature (at the lowest order) in a particular point of the 
space is 
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^m = U d\^-h rr . (4) 

J z J nuUgeodesic u 'l 

The integral in eq. ([4]) is taken over a path of null geodesic which leaves the 
current spacetime point heading off in the direction defined by the unit vector 
ft and going back to the surface of last scattering of the CMBR. 

Here A is a particular choice of the affine parameter along the null geodesic. 
By using conformal coordinates, we have for the metric perturbation 

5g a b = R 2 (rj)h ab , (5) 

and r in eq. (@]) is a radial spatial coordinate which goes outwards from the 
current spacetime point. The effect of a long wavelenght GW is to shift the 
temperature distribution of CMBR from perfect isotropy. Because the flucta- 
tions are very small (< AT/T Q >< 5 * 10~ 5 (221 ES]), the perturbations caused 
by the relic SBGWs cannot be too large. 

The WMAP results give rather tigh constraints on the spectrum of the 
SBGWs at very long wavelenghts. In [24} [25] we find a constraint on £l gw (f) 
which derives from the COBE observational limits, given by 

tt g w(f)hi Q0 < 7 * KT n (^) 2 for H < f < 30H . (6) 

Now the same constraint will be obtained from the WMAP data j2Q|. Be- 
cause of its specific polarization properties, relic SBGWs should generate par- 
ticular polarization pattern of CMBR anisotropies, but the detection of CMBR 
polarizations is not fulfiled today [31]. Thus an indirect method will be used. 
We know that relic GWs have very long wavelenghts of Hubble radius size, so 
the CMBR power spectrum from them should manifest fast decrease at smaller 
scales (hight multipole moments) . But we also know that scalar modes produce 
a rich CMBR power spectrum at large multipole moments (series of acoustic 
peaks, ref. |22l 23J). Then the properties of tensor modes of cosmological per- 
turbations of spacetime metric can be extract from observational data using 
angular CMBR power spectrum combined with large scale structure of the Uni- 
verse. One can see (fig. 1 ) that in the range 2 < / < 30 (the same used 
in [25], but with the old COBE data |32j) scalar and tensor contributions are 
comparable. From [22], 123) . the WMAP data give for the tensor/scala ratio r 
the constraint r < 0.9. (r < 0.5 in the COBE data, ref. [32]); Novosyadly and 
Apunevych obtained £l SC aiar{Ho) < 2.7 * 10 -9 |31j . Thus, if one remembers 
that, at order of Hubble radius, the tensorial spectral index is — 4 < n t < —2 
[20] . it results 

*W/)frioo < 1-6 * 10~ 9 (y) 2 for H < f < 30H , (7) 

which is greater than the COBE data result of eq. (jSJ. 
We emphasize that the limit of eq. (7| is not a limit for any GWs, but only 
for relic ones of cosmological origin, which were present at the era of the CMBR 
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Figure 1: The tensor to scalar ratio 



decoupling. Also, the same limit only applies over very long wavelenghts (i.e. 
very low frequencies) and it is far below frequencies of the Virgo - LIGO band. 

The primordial production of the relic SBGWs has been analyzed in [25l 
[26l [27l [33] (note: a generalization for f(R) theories of gravity will be given 
in the next section following [2]), where it has been shown that in the range 
10~ 15 Hz < f < 10 10 Hz the spectrum is flat and proportional to the ratio 

Pds „ 10 -i2 (g) 



PPlanck 

WMAP observations put strongly severe restrictions on the spectrum, as 
we discussed above. In fig. 2 the spectrum il gw is mapped: the amplitude 
(determined by the ratio ^ ds - ) has been chosen to be as large as possible, 
consistent with the WMAP constraint ((7j). 
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Figure 2: The spectrum of relic gravitons in inflationary models is flat over a 
wide range of frequencies. The horizontal axis is log 10 of frequency, in Hz. The 
vertical axis is log 10 O sr . The inflationary spectrum rises quickly at low frequen- 
cies (wave which rentered in the Hubble sphere after the Universe became mat- 
ter dominated) and falls off above the (appropriately redshifted) frequency scale 
fmax associated with the fastest characteristic time of the phase transition at 
the end of inflation. The amplitude of the flat region depends only on the energy 
density during the inflationary stage; we have chosen the largest amplitude con- 
sistent with the WMAP constrain discussed earlier: n gr (f )h'f Q0 < 1.6* 10~ 9 at 
10~ 18 Hz. This means that at Virgo and Lisa frequencies, Q gr (f)h 2 00 < 9*10 -13 

Nevertheless, because the spectrum falls off oc f~ 2 at low frequencies [25} 
EH E3 [33] , this means that today, at Virgo and LISA frequencies, indicated in 



%■ 2, 



*W/)fc?oo < 9*10 



(9) 



while using the COBE data it was 



^gw(f)h 2 W0 < 8 * l(T 14 (refs. [25l [32]). 



C 



It is interesting to calculate the correspondent strain at « 100Hz, where 
interferometers like Virgo and LIGO have a maximum in sensitivity. The well 
known equation for the characteristic amplitude (2Ql EH EH E2] can be used: 

h c (f) ~ 1.26 * 10- 18 (^) v // l ? 00 ^ tu (/), (10) 

obtaining 

/^(lo&ff*) < 1.7* icr 26 . (ii) 

Then, because for ground-based interferometers a sensitivity of the order of 
10~ 22 is expected at « lOOiTz, four order of magnitude have to be gained in 
the signal to noise ratio [9] SSI El S3 S3 SI S3 HE] ■ Let us analyze smaller 
frequencies too. The sensitivity of the Virgo interferometer is of the order of 
10~ 21 at w 10Hz OSQl and in that case it is 

ft c (10ffjz) < 1.7* 1(T 25 . (12) 

For a better understanding of the difficulties on the detection of the SBGWs 
a lower bound for the integration time of a potential detection with advanced 
LIGO is released. For a cross-correlation between two interferometers the signal 
to noise ratio (SNR) increases as 

(13) 

where is the one-sided power spectral density of the i detector [34] 

and 7(/) the well known overlap-reduction function [34] [35]. Assuming two 
coincident coaligned detectors (j(f) — 1) with a noise of order 10~ 48 /Hz (i.e. 
a typical value for the advanced LIGO sensitivity [36]) one gets (SNR) ~ 1 for 
~ 3*10 5 years using our result VL gw (f)h\ m ~ 9*10~ 13 while it is (SNR) ~ 1 for 
~ 3 * 10 7 years using previous COBE result VL gw (f)h\ m ~ 8 * 10~ 14 . Since the 
overlap reduction function degrades the SNR, these results can be considered 
a solid upper limit for the advanced LIGO configuration for the two different 
values of the spectrum. 

The sensitivity of the LISA interferometer will be of the order of 10 -22 at 
10 _ Hz [37] and in that case it is 

h c (10~ 3 Hz) < 1.7 * 1CT 21 . (14) 

Then a stochastic background of relic gravitational waves could be in prin- 
ciple detected by the LISA interferometer. We also hope in a further growth in 
the sensitivity of advanced projects. 

We emphasize that the assumption that all the tensorial perturbation in the 
Universe are due to a SBGWs is quit strong, but our results JSJ), (jlip . fT2|) and 
(SI]) can be considered like upper bounds. 

Reasuming, in this section the SBGWs has been analyzed with the auxilium 
of the WMAP data, while previous works in literature, used the old COBE 
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data, seeing that the predicted signal for these relic GWs is very weak. From 
our analysis it resulted that the WMAP bound on the energy spectrum and on 
the characteristic amplitude of the SBGWs are greater than the COBE ones, but 
they are also far below frequencies of the earth-based antennas band. In fact the 
integration time of a potential detection with advanced interferometers is very 
long, thus, for a possible detection we have to hope in a further growth in the 
sensitivity of advanced ground based projects and in the LISA interferometer. 

3 Tuning the stochastic background of gravita- 
tional waves with f(R) theories of gravity 

GWs are the perturbations h a /3 of the metric g a /3 which transform as three- 
tensors. Following ED [38], the GW-equations in the TT gauge are 



Dhj = 0, (15) 

where □ = {—g)~ 1 / 2 d a {—g) 1 / 2 g a Pdp is the usual d'Alembert operator and these 
equations are derived from the Einstein field equations deduced from the Hilbert 
Lagrangian density L = R [6l El] ■ Clearly, matter perturbations do not appear 
in (JH) since scalar and tensor perturbations do not couple with tensor pertur- 
bations in Einstein equations. The Latin indices run from 1 to 3, the Greek 
ones from to 3. Our task is now to derive the analog of eqs. lfT5|) assuming a 
generic theory of gravity given by the action 



1 

2^ 



d 4 x^f(R), (16) 



where, for a sake of simplicity, we have discarded matter contributions. A 
conformal analysis will help in this goal. In fact, assuming the conformal trans- 
formation 

g a p = e 2 *g Q /3 (17) 
where the conformal rescaling 



2* 



f'(R) (18) 



has been chosen being the prime the derivative with respect to the Ricci 
curvature scalar and $ the "conformal scalar field", we obtain the conformally 
equivalent Hilbert-Einstein action 

where i(<f>, $ ;Q ) is the conformal scalar field contribution derived from 

R a p = Rap + 2(* ;Q $;/3 - g a ^,8^ s - ^g af 3^ s ,8) (20) 
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and 

R=e~ 2 * + (R-6[]<f>-6<f> ;S & s ). (21) 

In any case, as we will see, the $ ;Q ,)-term does not affect the GWs-tensor 
equations so it will not be considered any longer (note: a scalar component in 
GWs is often considered \T7\ \39\ 140) 141] , but here we are taking into account 
only the genuine tensor part of stochastic background) . 

Starting from the action fl9|) and deriving the Einstein-like conformal equa- 
tions, the GWs equations are 

niv> = 0, (22) 

expressed in the conformal metric g a /3- Since no scalar perturbation couples 
to the tensor part of gravitational waves it is 

^'=ff% = e- 2 V j e%=4 (23) 
which means that h\ is a conformal invariant. 

As a consequence, the plane wave amplitude h\ = h{t)e\ exp(ifcjX J ), where 
el is the polarization tensor, are the same in both metrics. In any case the 
d'Alembert operator transforms as 

□ = e- 2 ' J> (a + 2fr a d., a ) (24) 

and this means that the background is changing while the tensor wave am- 
plitude is fixed. 

In order to study the cosmological stochastic background, the operator (f24|) 
can be specified for a Friedman- Robertson- Walker metric [Ml El], and the equa- 
tion 11221) becomes 

h+(3H + 2^)^ + k 2 a- 2 h = 0, (25) 
at at 

being □ = ^ + 3ff Jj, a(t) the scale factor and k the wave number. It has 
to be emphasized that equation ([251 applies to any f(R) theory whose conformal 
transformation can be defined as e 2 * = f{R). The solution, i.e. the GW 
amplitude, depends on the specific cosmological background (i.e. a(t)) and the 
secific theory of gravity (i.e. &(t)). For example, assuming power law behaviors 
for a(t) and = \ In f(R(t)), that is 

m)=f'(R)=ti(±-r, a(t) = a (±-r, (26) 
it is easy to show that general relativity is recovered for m = while 

m 2 + m - 2 

n = ^ — (27) 

m + 2 v ' 

is the relation between the parameters for a generic f(R) = foR" where 
s = 1— 'y with s \= 1 [42] . Equation (|25|) can be recast in the form 

h + (3n + m)r x ^ + k 2 a (j) 2n h = 0, (28) 
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whose general solution is 



h(t) = {jy^CMx) + C 2 J- a (x)}. (29) 

J a 's are Bessel functions and 

1 — in — m 1 — 3n — m kt l ~ n 
a = 2(n-l) ' P = 2 > X = — (30) 

while to , Ci, and C2 are constans related to the specific values of n and m. 

The time units are in terms of the Hubble radius H -1 ; n= 1/2 is a radiation- 
like evolution; n = 2/3 is a dust-like evolution, n = 2 labels power-law inflation- 
ary phases and n = —5 is a pole-like inflation. From eq. (|27|) . a singular case 
is for m = —2 and s = 2. It is clear that the conformally invariant plane- wave 
amplitude evolution of the tensor GW strictly depends on the background. 

Let us now take into account the iusse of the production of relic GWs con- 
tributing to the stochastic background. Several mechanisms can be considered 
as cosmological populations of astrophysical sources [43], vacuum fluctuations, 
phase transitions [24] and so on. In principle, we could seek for contributions 
due to every high-energy physical process in the early phases of the Universe 
evolution. 

It is important to distinguish processes coming from transitions like inflation, 
where the Hubble flow emerges in the radiation dominated phase and process, 
like the early star formation rates, where the production takes place during the 
dust dominated era. In the first case, stochastic GWs background is strictly re- 
lated to the cosmological model. This is the case we are considering here which 
is, furthermore, also connected to the specific theory of gravity. In particular, 
one can assume that the main contribution to the stochastic background comes 
from the amplification of vacuum fluctuactons at the transition between an in- 
flationary phase and the radiation-dominated era. However, in any inflationary 
model, we can assume that the relic GWs generated as zero-point fluctuations 
during the inflation undergoes adiabatically damped oscillations (~ l/o) until 
they reach the Hubble radius H~ 1 . This is the particle horizon for the growth of 
perturbations. On the other hand, any other previous fluctuation is smoothed 
away by the inflationary expansion. The GWs freeze out for a/k~3> H^ 1 and re- 
enter the H^ 1 radius after the reheating in the Friedman era [2l" 1 l25l l26l Wf\ [33] . 
The re-enter in the radiation dominated or in the dust-dominated era depends 
on the scale of the GW. After the re-enter, GWs can be detected by their Sachs- 
Wolfe effect on the temperature anisotropy at the decoupling [30]. When $ 
acts as the inflaton |21l [44] we have ^? <C H during the inflation. Considering 
also the conformal time dr/ = dt/a, eq. lf25|) reads 

h" + 2—ti + k 2 h = 0, (31) 

7 

where 7 = ae*and derivation is with respect to 77. Inflation means that 
a(t) = clq exp(Ht) and then r\ = J r dt/a = l/(aH) and ^ = — ~. The exact 
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solution of Il31|) is 

h(r]) = k^ 3/2 ^/2/k[C 1 (sm kr/ - cos kr]) + C 2 (sm kr] + cos kr/)}. (32) 

Inside the 1/H radius it is kr] ^> 1. Furthermore considering the absence of 
gravitons in the initial vacuum state, we have only negative-frequency modes 
and then the adiabatic behavior is 

h = k l/2 ^/2j^— Cexp(-ikrj). (33) 
aH 

At the first horizon crossing (aH = k), the averaged amplitude Ah = 
(k/2n) 3 / 2 \h\ of the perturbations is 

A h = (34) 

when the scale a/k grows larger than the Hubble radius 1/H, the growing 
mode of evolution is constant, that it is frozen. This situation corresponds to the 
limit — kr] <C 1 in equation ((32)) . Since $ acts as the inflaton field, it is $ ~ 
at re-enter (after the end of inflation). Then the amplitude Ah of the wave 
is preserved until the second horizon crossing after which it can be observed, 
in principle, as an anisotropy perturbation of the CBR. It can be shown that 
^rjr- < Ah is an upper limit to Ah since other effects can contribute to the 
background anisotropy [45]. From this consideration, it is clear that the only 
relevant quantity is the initial amplitude G in equation ([33| which is conserved 
until the re-enter. Such an amplitude directly depends on the fundamental 
mechanism generating perturbations. Inflation gives rise to processes capable 
of producing perturbations as zero-point energy fluctations. Such a mechanism 
depends on the adopted theory of gravitation and then could constitute a 
further constraint to select a suitable f(R)-theory. Considering a single graviton 
in the form of a monocromatic wave, its zero-point amplitude is derived through 
the commutation relations 

[h(t,x),ir h (t,y)\ = iS 3 (x-y) (35) 

calculated at a fixed time t, where the amplitude h is the field and irh is the 
conjugate momentum operator. Writing the lagrangian for h 

L = \^/tg al3 h, a h,p (36) 
in the conformal FRW metric g a f3 (h si conformally invariant), we obtain 



Equation (|35|) becomes 



r, / % dh , s , S 3 (x — y) ,„„. 
[h(t,x),-(y,y)]=i ^J> (38) 
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and the fields h and 4r can be expanded in terms of creation and annihilation 
operators 

h(t,x) = J d'k\h{t)e- ik > + A-(i)e' fa ] (39) 

The commutation relations in conformal time are then 

[W* - = ( 4l ) 

Insering ([33]) and ([3"4"]l. we obtain C = \/27r 2 i/e~* where i? and $ are 
calculated at the first horizon crossing and then 

A h = ^He-\ (42) 

which means that the amplitude of GWs produced during inflation directly 
depends on the given f(R) theory being <I> = \ In f'(R)- Explicitly, it is 

A h = - H . (43) 

This result deserves some discussion and can be read in two ways. From one 
side the amplitude of relic GWs produced during inflation depends on the given 
theory of gravity that, if different from general relativity, gives extra degrees 
of freedom which assume the role of inflaton field in the cosmological dynamics 
[44] . On the other hand, the Sachs- Wolfe effect related to the CMBR temperatue 
anisotropy could constitute a powerful tool to test the true theory of gravity at 
early epochs, i.e. at very high redshift. This probe, related with data a medium 
|46| and low redshift |47j . could strongly contribute 

1. to reconstruct cosmological dynamics at every scale; 

2. to further test general relativity or to rule out it against alternative theo- 
ries; 

3. to give constrains on the SBGWs, if f(R) theories ares independently 
probed at other scales. 

Reasuming, in this section it has been shown that amplitudes of tensor GWs are 
conformally invariant and their evolution depends ond the cosmological SBGWs. 
Such a background is tuned by conformal scalar field which is not present in 
the standard general relativity. Assuming that primordial vacuum fluctuations 
produce a SBGWs, beside scalar perturbations, kinematical distorsions and so 
on, the initial amplitude of these ones is function of the f(R)-theory of gravity 
and then the SBGWs can be, in a certain sense, "tuned" by the theory. Vice 
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versa, data coming fro the Sachs-Wolfe effect could contribute to select a suit- 
able f(R)-theory which can be consistently matched with other observations. 
However, further and accurate studies are needed in order to test the relation 
between Sachs- Wolfe effect and f(R) gravity. This goal could be achieved in the 
next future through the forthcoming space (LISA) and ground based (Virgo, 
LIGO) interferometers. 



4 Conclusions 

The SBGWs has been analyzed with the auxilium of the WMAP data while, 
in general, in previous works in the literature about the SBGWs, old COBE 
data were used. After this, it has been shown how the SBGWs and f(R) gravity 
can be related, showing, vice versa, that a revealed SBGWs could be a powerly 
probe for a given theory of gravity. In this way, it has also been shown that the 
conformal treatment of SBGWs can be used to parametrize in a natural way 
f(R) theories. 
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